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Abstract 

The  response  of  a unity-feedback  system  with  a delay  element  in  the  forward  path  exhibits 
a periodic  component  that  can  be  approximated  by  truncating  its  harmonic  expansion. 
Rational  approximants  of  the  transfer  function  e~Ts  of  such  element  can  simply  be  ob- 
tained from  this  closed-loop  approximation.  A unifying  approach  to  recent  methods  based 
on  this  criterion  [2,  3]  is  presented,  which  allows  us  to  point  out  their  respective  features. 
The  standard  Pade  technique  and  a heuristic  method  described  in  [5]  are  also  considered. 

1 Introduction  and  problem  statement 

In  modelling  dynamic  systems  for  control  purposes,  it  is  often  necessary  to  account  for 
time  delays  due,  e.g.,  to  transport  phenomena  or  distributed-parameter  components. 

The  response  of  an  ideal  delay  element  (delayor)  to  an  input  u(t),  identically  equal  to 
0 for  t < 0,  is  y(t)  = u(t-T ),  T > 0,  where  T indicates  the  time  delay.  By  denoting  with 
U(s)  the  Laplace  transform  of  u(t),  the  Laplace  transform  of  y(t)  is  Y(s)  = e~TsU(s). 
Therefore  the  transfer  function  of  the  delayor  is  the  transcendental  function  e~Ts. 

The  problem  of  approximating  e~Ts  by  means  of  a rat  ional  function  has  a long  history 
(see,  e.g.,  [1])  but  is  still  important  from  both  the  computational  and  the  conceptual 
point  of  view;  a few  recent  contributions  on  the  subject  are  quoted  in  [2].  In  many 
practical  applications  the  physical  realizability  and  the  stability  of  the  approximant 
limit  the  choice  of  the  approximant  to  proper  rational  functions  with  real  coefficients 
and  a Hurwitz  denominator.  These  requirements  are  satisfied  by  Blaschke  products,  i.e., 
functions  of  the  form: 

B ^ = nr=i(s  + ai)'  Ke[ai]>0.  (1.1) 
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This  has  the  desirable  property  that  \B(ju)\  = |e_jT“|  = 1,  Vw,  and  arg[B(jw)]  is 
monotonically  decreasing  with  oj  like  arg[e~jT“]  = — Tw.  On  the  other  hand,  the  step 
response  of  a system  with  transfer  function  B(s)  starts  from  +1  or  -1,  whereas  the  step 
response  of  an  ideal  delayor  obviously  starts  from  0. 

The  most  widely  adopted  method  to  form  a rational  approximant  of  a delay  element  is 
based  on  the  Pade  technique  which  does  not  always  guarantee  stability  (even  if  biproper 
Pade  models  are  necessarily  stable).  Since  such  a technique  leads  to  the  retention  of  the 
first  Maclaurin  expansion  coefficients  of  e~Ts,  the  resulting  approximation  is  the  best  in 
the  neighbourhood  of  uj  = 0.  In  different  frequency  bands,  other  types  of  models  may  be 
preferred. 

In  [3]  a unity-feedback  system  whose  forward  path  consists  of  a delayor  is  analysed. 

In  the  case  of  negative  feedback,  the  unit  step  response  is  a piecewise  constant  function 
taking  on  the  value  0 for  2 kT  < t < (2k  + 1)T  and  the  value  1 for  (2k  + 1 )T  < t < 
(2k  + 2)T,  k > 0,  which  can  be  decomposed  into  a step  of  amplitude  |,  and  a square 
wave  of  amplitude  \ starting  from  — | at  t = 0. 

In  the  case  of  positive  feedback,  similar  considerations  allow  us  to  decompose  the  unit 
step  response  into  a linear  ramp  of  slope  a step  of  amplitude  — and  a saw-tooth 
wave  that  linearly  decreases  from  \ to  in  every  period  from  kT  to  (k  + 1)T. 

In  both  cases,  the  periodic  component  can  easily  be  expressed  as  a series  of  har- 
monic terms  (for  t > 0).  It  is  therefore  natural  to  approximate  the  step  response  of  the 
unity-feedback  system  by  retaining  the  non-periodic  component  together  with  a suitable 
number  of  the  first  harmonics  of  the  periodic  component. 

A rational  approximation  Wa(s)  of  the  transcendental  transfer  function  W(s)  of  the 
above-mentioned  feedback  system  is  obtained  by  dividing  the  Laplace  transform  of  the 
approximate  step  response  by  the  Laplace  transform  j of  the  step  input.  The  rational 
approximant  Ga(s)  of  the  delayor  transfer  function  is  then  determined  as 


Ga(s)  = 


Wa(s) 
lTW^s)  ’ 


(1.2) 


where  the  minus  sign  applies  to  the  case  of  negative  feedback  and  the  plus  sign  to  that 
of  positive  feedback.  It  turns  out  [3]  that  Ga(s)  is  a stable  biproper  rational  function 
having  the  form  of  a Blaschke  product;  precisely,  negative  feedback  supplies  even-order 
approximants  and  positive  feedback  produces  odd-order  approximants. 

Obviously,  the  same  result  could  be  achieved  by  referring  to  different  inputs  (even 
an  impulse),  but  the  choice  of  the  unit  step  is  particularly  convenient.  According  to 
the  terminology  suggested  in  [4],  the  rationale  of  such  a procedure  consists  in  retaining 
the  “input  component”  (and  the  “resonant  component” , if  any)  and  in  truncating  the 
periodic  “system  component”  of  the  response. 

In  [2]  a feedback  structure  is  used  as  well,  but  another  approximation  criterion  is 
adopted,  which  leads  to  different  models  depending  on  the  chosen  input.  In  particular, 
the  family  of  inputs  considered  in  [2]  is  {u(t)  = tm  ,m  £ N ,t  > 0},  and  the  procedure 
exploits  several  properties  of  Bernoulli  numbers  and  polynomials. 

In  the  following,  the  above  approaches  are  presented  in  a unified  form  which  allows 
us  to  point  out  their  respective  features  and  to  derive  the  related  approximants  in  an 
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easier  way.  Finally,  criteria  are  given  to  choose  the  approximation  that  is  most  suited  to 
the  application  at  hand,  also  taking  into  account  the  standard  Pade  approximation  and 
a further  approximation  presented  in  [5]. 

2 Derivation  of  the  approximant 

For  the  sake  of  simplicity,  we  shall  almost  exclusively  refer  to  the  case  of  negative  feed- 
back; only  a brief  mention  will  be  made  of  the  case  of  positive  feedback. 


2.1  Negative  feedback 

The  transfer  function  W ( s ) of  the  negative  feedback  system  with  forward-path  transfer 
function  G(s ) = e~Ts  is 

<21> 

whose  singularities  (poles)  are  the  roots  of  eTs  = — 1 , i.e. 

s-  ±jpk  :=  ±j(2k  - 1)^,  k € Z+  . 

IF(s)  can  also  be  interpreted  as  the  Laplace  transform  of  the  sequence  of  positive  and 
negative  impulses  forming  the  derivative  of  the  step  response  described  in  the  intro- 
duction. Therefore,  it  is  the  sum  of  a constant  equal  to  \ (corresponding  to  the  step 
component  in  the  just-mentioned  step  response)  and  a series  of  “harmonic”  terms  asso- 
ciated with  the  above  poles: 

W(s)  = l + f]  |-^r-  + — — ] , 

2 s-jpk  s + jpk 

where  the  bar  denotes  conjugate  and,  using  the  standard  formula  for  the  residues, 


Tk  — lim  (s-jpk)W(s) 

S->JPk 


It  follows  that 


i 9 00 


fr{s^  + (2k- 1)2  £ 


In  order  to  compare  the  results  in  [2]  and  [3],  let  us  consider  a canonical  input  of  the 
form 

ui(t)  = jr  — i t > 0,  (2-3) 


whose  Laplace  transform  is 


Ui(s)  = Ti. 


(In  [3]  only  the  case  of  z — 1 is  considered,  whereas  the  inputs  used  in  [2]  differ  from 
(2.3)  by  a scaling  factor  which  is  irrelevant  for  the  following  considerations.) 


Remarks  on  delay  approximations  based  on  feedback 


415 


On  the  basis  of  (2.2)  the  Laplace  transform  of  the  (forced)  response  to  (2.3), 


yi(5.)  = _ W(s) , 


can  be  rewritten  as 


where  for  i even, 


and  for  i odd, 


i— 1 


Y(s)  = -Y  chsh  + V Uki  + ^ kiS 


h= 0 


Oiki  — 0 5 


Phi  = ("I)5 


Oiki  = (-1)  2 


<_2 

T [(2/c  — 1)tt 

(i-i) 


Pki  = 0 . 


T L(2fc-  1)tt 

Therefore,  IT (s)  can  also  be  presented  in  the  alternative  form 

yes')  . .22.  4- 

W{a)  - *W 

(L'l.M 

h=0 


Ui(s) 


h , Qfc*s*  + 


(2.40 


(2.4n) 


(2.5) 


Each  term  of  the  series  in  (2.5)  is  given  by  the  sum  of  a polynomial  of  degree  i — 1 
(quotient  of  the  division  of  its  numerator  by  its  denominator)  and  a strictly  proper 
rational  function  (whose  numerator  is  the  remainder  of  the  division).  Therefore,  (2.5) 
becomes 

2—1  OO  ( 2—1  i X ^ 

fc=l  U=o 


h= 0 

which  can  be  rewritten  as 
2 — 1 


s2  + (2fc-l)2f| 


(2.6) 


W(s)  = E Ch  + E^K+E 


Ifki  + dkiS 


h= 0 


fc=l 


^s2  + (2fc-l)2fi 


By  comparing  (2.7)  with  (2.2),  one  finds  that 

OO 

CO  + dki, 0 

k= 1 

CO 


2 ’ 


c-h  T ^ ^ dkith  — d , h > 0, 


(2.7) 

(2.8) 

(2.9) 


fc=i 


7ki 


0,  Vfc,f, 


^ , V k,  i . 

The  procedure  suggested  in  [3]  could  alternatively  be  presented  with  reference  to 
expression  (2.7)  where  coefficients  related  to  the  specific  input  appear.  Precisely,  the 
approximant  W„(s)  is  obtained  in  this  case  by  adding  to  the  exact  value  \ of  the  first 
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sum  (cf.  (2.8)  and  (2.9))  the  first  K (harmonic)  terms  of  the  second  summation 

w (s)  = - - - Y - - 

which  is  independent  of  the  input  u,(f). 

The  procedure  suggested  in  [2]  refers  instead  to  expressions  (2.5)  or  (2.6),  and  the 
approximation  consists  in  truncating  the  summation  over  k,  where  each  addendum  is 
formed  by  a polynomial  and  a strictly  proper  harmonic  term.  Therefore  the  resulting 
Wa(s ) is 

h= 0 fc=l  S + T* 

which  does  depend  on  i and  it  is  not  proper  because  the  part  added  to  the  harmonic 
terms  does  not  reduce  to  the  constant  as  is  instead  the  case  in  IT(s).  Nevertheless, 
the  approximant  Ga(s)  — W„(s)/(1  — H/C[(.s))  of  e~Ts  turns  out  to  be  biproper. 

As  concerns  the  computation  of  the  above  approximants,  the  suggested  approach 
seems  to  be  preferable  to  that  adopted  in  [2]  because 
(i)  coefficients  Ch,  which  correspond  to  the  first  i Maelaurin  expansion  coefficients  of 

^ h\ 


can  be  easily  be  evaluated  using  the  classic  Pade  procedure,  and 
(ii)  formulae  (2.4*)  and  (2.4ii)  immediately  supply  coefficients  ctki,f3ki- 

2.2  Positive  feedback 

Considerations  analogous  with  those  of  Section  2.1  lead  to  the  following  transfer  function 
in  the  case  of  positive  feedback 


119 

"'w-K-3 (2-n) 

so  that  Yi(s)  = VP(s)t/j(s)  can  be  separated  into  a (harmonic)  series  associated  with  the 
imaginary  conjugate  poles  of  IT(s)  and  a strictly  proper  fraction  with  denominator  .s,+1 . 
Using  the  terminology  in  [4],  the  mentioned  series  corresponds  to  the  “system  compon- 
ent” of  the  forced  response  and  the  fraction  corresponds  to  its  “interaction  component” 
because  the  poles  of  the  latter  are  common  to  !U(s)  and  E/,(s)  (no  “input  component” 
is  present  in  this  case  since  Ui(s)  does  not  exhibit  poles  different  from  those  of  W(s)). 

As  shown  in  [3],  the  truncation  of  the  series  in  (2.2)  results  in  even-order  biproper 
approximants  Ga(s),  whereas  the  truncation  of  the  series  in  (2.11)  results  in  odd-order 
biproper  approximants  Ga{s). 

Instead,  as  shown  in  [2],  truncating  the  series  in  (2.5)  leads  to  odd-order  approxim- 
ants, whereas  truncating  the  analogous  series  corresponding  to  positive  feedback  leads 
to  even-order  approximants. 
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2.3  Stability  and  approximation  error 


It  has  been  proved  [3]  that  the  even-order  rational  approximations  Ga(s)  of  e~Ts  obtained 
from  (2.1),  as  well  as  the  odd-order  ones  obtained  by  truncating  (2.11),  are  stable. 
Instead,  as  explicitly  stated  in  [2]  for  inputs  tm,  m > 2 (i.e.,  using  the  previous  notation, 
Ui(t)  with  i > 3)  the  “alternating  sign  of  the  Bernoulli  numbers  makes  the  approximation 
in  general  unstable  [...].  Hence,  from  a practical  point  of  view,  any  improvement  with 
respect  to  the  approximants  obtained  in  [3]  is  to  be  found  with  p — 1”,  i.e.,  i = 2, 3. 


The  approximation  accuracy  can  be  evaluated  by  referring,  e.g.,  to  the  “closed-loop 
error” 

E(s)  :=  W(s)  - Wa(s ) . 

From  (2.1)  we  get 

c\  °° 

E(s)  = E1(s):=-~  £ 

k=K+ 1 

whereas  from  (2.10)  we  have 


2 ^ s 
T . 4^  . s2  + pi  ' 


i— 1 oo 

E(s)  = E2(s):=Y1  £ 


dki,h.Sh  + E\(s)  . 


h= 0 k=K+ 1 


Since  E(s)  is  a complex  quantity,  | JE7a(s) [ may  well  be  smaller  than  |£u(s)|  for  certain 
values  of  $ (or  jw ) . 


3 Alternative  approximants 


As  already  pointed  out,  the  procedure  suggested  in  [2]  leads  to  approximants  that  depend 
on  the  chosen  canonical  input.  To  improve  the  approximation  within  suitable  frequency 
bands  not  centred  at  the  origin,  it  is  reasonable  to  resort  to  non-canonical  inputs  whose 
spectrum  has  larger  amplitude  there.  A simple  choice  corresponds,  e.g.,  to 


U(s) 


1 + 2£~—  + —5- 


U)ri 


in  which  ton  is  at  the  centre  of  the  band  and  £ is  suitably  small. 


The  choice  of  the  form  of  the  input  (as  well  as  the  order  of  the  canonical  input)  is 
somewhat  arbitrary  and  is  influenced,  in  practice,  by  empiric  considerations.  Therefore, 
it  makes  sense  to  compare  the  results  of  the  above  procedures  with  those  obtained  in 
[5]  using  a heuristic  procedure  based  on  the  direct  approximation  of  the  phase  Bode 
diagram  of  by  means  of  a Blaschke  product  Bn(ju> ) of  order  n.  For  n odd,  the 

first  factor  of  Bn(s)  has  the  form 


1 — TS 


and  the  others  have  the  form 


1-2&— +* 


Gi(s) 


UJtj 


W‘ 


„ s s 

1 + 2 1 Y 


UJn 


(Ol 


1 > ^ > 0 , (Oni  ^ 0 > 


(3.1) 


whereas  for  n even  all  factors  have  form  (3.1). 

All  the  considered  techniques  produce  unit-magnitude  all-pass  frequency  responses  so 
that  the  approximation  they  afford  can  be  judged  with  reference  to  the  phase  deviation 
A(joj)  from  —Tu)  only.  As  uj  — > oo,  A (jui)  — ► oo  in  all  cases.  Therefore,  reasonable  criteria 
for  choosing  the  method  most  suited  to  the  specific  application  are:  (i)  the  bandwidth 
Be  where  |A(jw)|  is  less  than  a specified  value  e,  or  (ii)  the  maximum  A b of  |A(jw)|  in 
a prescribed  band  B. 

By  way  of  example,  Fig.  1 shows  A (ju)  vs  ui  for  the  4-th  order  all-pass  approximants 
of  e~iu,  (T  = 1)  obtained  according  to  (2.1)  with  K — 2 (curve  a),  to  the  procedure 
suggested  in  [2]  for  Us(t)  = t2  (curve  b),  to  the  standard  Pade  procedure  (curve  c),  and  to 
the  heuristic  method  in  [5]  (curve  d).  For  instance,  with  reference  to  criterion  (i)  above, 
the  Pade  approximant  is  best  for  e very  small,  the  method  suggested  in  [2]  is  optimal 
for  e ~ 10°,  and  the  heuristic  method  and  the  method  suggested  in  [3]  are  preferable  for 
£ > 45°. 

Analogous  results  are  obtainable  for  approximants  of  different  order. 


Fig.  1.  Phase  deviations  A (ju)  for  the  considered  4-th  order  approximants. 


4 Conclusions 

The  approximation  procedure  presented  in  [2]  and  [3]  have  been  embedded  in  a unified 
frame  which  points  out  well  their  respective  features  and  allows  us  to  determine  the 
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parameters  of  the  approximants  in  an  easier  way.  Criteria  have  been  provided  for  choosing 

the  approximation  method  that  is  most  suited  to  the  specific  application. 

Bibliography 

1.  0.  Perron,  Die  Lehre  von  den  Kettenbriichen.  Stuttgart:  Teubner,  1913.  3rd  ed. 

1957.  In  German. 

2.  C.  Battle  and  A.  Miralles,  “On  the  approximation  of  delay  elements  by  feedback,” 
Automatical  vol.  36,  pp.  659-664,  2000. 

3.  A.  Beghi,  A.  Lepschy,  and  U.  Viaro,  “Approximating  delay  elements  by  feedback,” 

IEEE  Trans.  Circ.  Sys.  I,  vol.  44,  pp.  824-828,  1997. 

4.  P.  Dorato,  A.  Lepschy,  and  U.  Viaro,  “Some  comments  on  steady-state  and  asymp- 
totic responses,”  IEEE  Trans.  Education,  vol.  37,  pp.  264-268,  1994. 

5.  A.  Beghi,  A.  Lepschy,  and  U.  Viaro,  “On  the  simplification  of  the  mathematical 
model  of  a delay  element,”  in  E.  Kuljanic,  ed.,  Advanced  Manufacturing  Systems  and 
Technology:  Springer  Verlag,  1996,  pp.  617-624. 


